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P ! ABSTRACT 

^ ' Wc present high-speed optical photometry of the soft gamma repeater 

^ I SGR 0501+4516, obtained with ULTRACAM on two consecutive nights approxi- 

mately 4 months after the source was discovered via its 7-ray bursts. We detect 
SGR 0501-1-4516 at a magnitude of i' = 24.4 ± 0.1. We present the first measure- 
ment of optical pulsations from an SGR, deriving a period of 5.7622 ± 0.0003s, in 
^ ' excellent agreement with the X-ray spin period of the neutron star. Wc compare the 

morphologies of the optical pulse profile with the X-ray and infrared pulse profiles; 
we find that the optical, infrared and harder X-rays share similar double-peaked mor- 
phologies, but the softer X-rays exhibit only a single-peaked morphology, indicative 
of a different origin. The optical pulsations appear to be in phase with the X-ray pul- 
'^O . sations and exhibit a root-mean-square pulsed fraction of 52 ± 7%, approximately a 

factor of two greater than in the X-rays. Our results find a natural explanation within 
the context of the magnetar model for SGRs. 
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1 INTRODUCTION ferred dipole magnetic field s trengths of B-10"-10^^G 

. (nowever, see iRea et al.|[2010t ). The loss of rotational en- 
Soft samma repeaters (SGRs) derive their name from . . m • . . 1 • .1 • . , ir 1 
, . , . , , , , ergy IS msuthcient to explain the persistent X-ray lumi- 
tne manner m which they are usually discovered - the .,11 1. .. . .1 

nosity and hence an alternative energy source is required, 

emission ot repeated short bursts ot hard X-rays/sott rr^, , , , , j 1 • ii, j r 

' ihe most-commonly accepted model is the decay ot the 

7-rays with durations of 0.1— 0.2 s and luminosities of , , . r- 1 1 mi • nit j , 11 

in38 in4i -1 t f -ti ar^T} -4- strong magnetic held, ihis so-called magnetar model can 

explain both the bursts/flares and the persistent emis- 



ergs . Less frequently, SGRs emit even more 

l41 -I n43 



energetic (10 —10 ergs ) bursts of longer duration . j. orrn r ■ 1 • r ii 

, , o y a gjQjj qj- ^i^g SGRs, the lormer via cracking 01 the neutron 

(1— 60s), known as intermediate flares, and rarer still, , . • 1 11 .1 1 -n- r- 1 1 1 .1 1 .. 

+ . n , /in44 in47 -1 onn cnn ^ A ■ u- u ^^^^ induced by the drifting field, and the latter 

giant fiares (10 —10 ergs , 300— oOOsJ, during which, 



via internal heating and mag netospheric currents induced 

tor a fleeting instant, the SGR outshines the entire i, ii, j -a^- 1^ ■ ^- n u — T~F. IFTT^ 

. ' | „, , mv. 1 by the driftmg/twistm g field (^Thompson fc DuncanI 11993. 

Galaxy, hor a review 01 SGKs, see I Woods & IhompsonI rn?r r^lnnnoK a li i- ii, ■ ^ ii 

ll 200d ) d [m — h ttil (I 2OO8I ) — I i hompson et al.l 120021 1 . Alternative theories to the magne- 

" ^' 7 ' — 7+ ;n n^T^ tar model d o exist, the most prominent being the 'fall-back 

As weU as periods of burstmg/flarmg, SGRs also ex- ^^^^^ ^yan pLradiis et al.lM . lGhatteriee et allB . 

hibit persistent X-ray emission, trom which it has been de- rr; ]|r,r.r>ih t • 1 1 j_i • ^ j_ v • • 

, , „ . , , . , , lAlpan 120011 1. In this model, the persistent X-ray emission 

duced that SGRs are all isolated neutron stars with rel- „ nr^r> 1 ii 11 1 j_ 1 a 1 ^ r> 1 

trom SGRs, and the closely-related Anomalous X-ray Pul- 
atively long spin periods m the range 2— 9 s, relatively / a -i^r> \ ■ 1 1 j_i j_- r j_ • 1 1 

. , . , ^ ,„_in -,^-12 -1 , , . n . sars (AXPs), is powered by the accretion ot material left 

rapid spm-down rates ot '^lO —10 ss and high m- r .1 1 . . r 1 .1 

over trom the supernova explosion that formed the neutron 



star. 
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Only 21 SGRs and AXPs are currently knowi{3, 7 of 
which are confirmed as SGRs (although the distinction be- 
twe en SGRs and A XPs is becoming increasingly blurred, 
e.g. iRea et al.|[2009l ). The discovery of each new source is 
therefore of considerable scientific interest. On 2008 August 
22, the Swift Burst Alert Telescope (BAT), observed mul- 
tiple, short (< 128 ms) bursts of soft (< 100 keV) 7-rays 
from a source close to the Galactic plane (|Barthelmv et all 
l2008h . The source's position, and the duration, energy, 
and repetition of the bursts, led to the classification of 
the object as a soft gamma repeater, SGR050H-4516 
(jBarthelmv et al]|2008l l: at the time, this was the first new 
SGR to be discovered for a decade. Subsequent observations 
of th e X-ray afterglow allowed the sp in period |Gogiis et al] 
l2008h and rate of period decrease ijWoods et al.l l2008h to 
be measured, further supporting the source's classifica- 
tion as an SGR. Detailed studies of SGR 05 1-1-45 16 in 
the X-rays have b een published by iRea et al] (|2009l ) and 
iGogiis et al] (I2010I). No radio counterpart of the SGR was 
detected ( Hessels et al.l [2OO8I ), but a faint, possibly vari- 
able, object coincident with the X-ray afterglow was dis- 
covered in the infrared jTanvir fc Varricattll2"oo3 ) and opti- 
cal (|Fatkhullin et al] I2OO8I ). In this paper, we present the 
results of follow-up, time-resolved optical observations of 
SGR 0501+4516, obtained with the aim of understanding 
the emission processes at work. 




Figure 1. i'-band finding chart for SGR 0501+4516. Left: 
Summed i'-band image of the field around SGR 0501+4516 on 
2009 January 1, with a total exposure time of 3247 s. The box 
shows the portion of the field that is plotted at a higher contrast 
on the right. The orientation of the image is marked on the up- 
per left-hand side. The pixel scale is 0.3 "/pixel, hence the field of 
view in this image is 60" X 60" . The vertical streaks are due to the 
frame transfer process. Right: Higher-contrast plot of a 16"x22" 
field around SGR 0501+4516, highlighting the detection of the 
pulsar in the i' band. The X-ray error circle lies within the seeing 
disc of the pulsar and is hence invisible. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed the proposed optical counterpart of 
SGR 0501+451 6 using the high-sp eed, triple-beam camera 
ULTRACAM (|Dhillon et al.l |2007^ on the 4.2 m William 
Herschel Telescope on La Palma. We obtained 10 152 g' 
and i' frames, each of 0.475 s exposure time (a run length 
of 1.4 hr), on 2008 December 31, and a further 6836 such 
frames (a run length of 1.0 hr) 26 hours later on 2009 
January 1. ULTRACAM was used in two- windowed mode, 
giving a dead time between each frame of 0.024 s, where 
each frame is time-stamped to a relative (i.e. frame-to- 
frame) accuracy of ~50/xs and to an absolute accuracy o f 
~lms using a dedicated GPS system (|Dhillon et al. 2007h. 
Obse rvations of the SDSS standard G 243-38 (|Smith et all 
|2002| ) were also obtained to flux calibrate the data. Both 
nights were photometric with seeing of 0.8" — 1.2" and 
first-quarter Moon. The sum of the i'-band frames obtained 
on 2009 January 1 is shown in Figure [T] 

The data were reduce d using the standard ULTRACAM 
data reduction software (|Dhillon et al.l l2007l l . All frames 
were first debiased and then flat-fielded, the latter using the 
median of twilight sky frames taken with the telescope spi- 
ralling. Adopting the same successful ap proach that we used 
in our studies of the AXPs 4U 0142+61 piiillon et al.lb'oosi ') 
and IE 1048.1-5937 (|DhiUon et al.ll200g ). we extracted light 
curves of SGR 0501+4516 using two different techniques: 



^ McGill SGR/AXP Online Catalogue: 

|http: / /www. physics. mcgill.ca/ ~pulsar/magnetar/main.html| 



2.1 Technique (i) 

First, we extracted a light curve of SGR 0501+4516 by shift- 
ing each frame to correct for image motion and then adding 
them into 10 evenly-spaced phase bins using the X-ra y 
epoch, spin period and spin-down rate of iRea et al.l ()2009h . 
This ephemeris is accurate enough to cover the epoch of 
our optical observations and is in go od agreement with the 
original X-ray ephemeris report ed bvllsrael et alj (|2008h . An 
optimal photometry algorithm (|Navloij | 19981 was then used 
to obtain sky-subtracted counts from SGR 0501+4516 and 
an i'~13 comparison star ~28" to the east of the SGR, 
the latter acting as the reference for the profile fits and 
transparency-variation correction (although this correction 
made only a negligible difference to the light curves). The 
position of SGR 0501+4516 relative to the comparison star 
was determined from a sum of all the images, and this offset 
was then held fixed during the reduction so as to avoid aper- 
ture centroiding problems. The sky level was determined 
from a clipped mean of the counts in an annulus surround- 
ing the target stars and subtracted from the object counts. 



2.2 Technique (ii) 

A light curve was also obtained using an independent 
method, where we omitted the phase binning step described 
in Section [2. II and simply performed optimal photometry on 
the 16 988 individual data frames, followed by a search for 
the X-ray period in a periodogram. 
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3 RESULTS 
3.1 Magnitudes 

Figure [1] shows that we have clearly detected an object 
in the i' band with magnitude i' = 24.4±0.1 at position 

RAj2ooo = 05'^01'"06.73%5j2ooo = 45°16'34.1"(±0.3"), in 
excellent agreem ent with the Cha ndra X-ray position of 
SGR 0501+4516 (|Gogiis et al.ll2010D . The source remained 
undetected in the u and g' bands to 3a limiting magnitudes 
of u' > 24.7 and g' > 26.9, as expected du e to the high hy - 
drogen column density {Nh ~ 10^^ cm~^. iRea et al.1 120091 ') 
at the low galactic latitude of the SGR. The i'-band mag- 
nitude of the source in our observations is fainter than the 
value s of 7 = 23.3 ±0.4 on 2008 August 25 (|Fatkhullin et al 
2008[) and i' = 23.5 ± 0.4 on 2008 September 2 (|Ofek et al 
20081 ). indicating that the SGR had declined in brightness by 
approximately one magnitude in the 4 months since ouburst. 



3.2 Optical light curves 

The first data reduction technique (Section 12. l|l resulted in 
the light curves shown as solid lines in Figures [^KilSb and[21:. 
The light curves clearly show a double-peaked pulse profile, 
with a broad peak around phase 0.15 and a narrower peak 
around phase 0.55 (where phase is arbitrarily defin ed as 
MJD 54701.0 in the X-ray ephemeris of |ReaelaL|[2Q0i) . The 
fact that the profile exhibits approximately the same mor- 
phology, phasing and pulsed fraction on both nights is strong 
evidence that the optical pulsations are real. 

An independent test of the reality of the optical pulsa- 
tions can be obtained from the second data reduction tech- 
nique (Section I2.2|l . Rather than folding the data on the 
X-ray ephemeris, we can instead determine the pulse pe- 
riod directly from our optical data using a periodogram. 
The non-sinusoidal nature of the pulsation profile led us to 
use the Phase Dis persion Minimization (PDM) technique 
llStellingwerj Il97^ ) to construct a periodogram, where the 
light curve is folded on a set of trial frequencies into 10 
phase bins and the normalised variance of the data in each 
binned light curve is then plotted as a function of frequency. 
The results are shown in Figure [3];. The deepest trough, 
and hence most likely period, in the periodogram occurs 
at 5.7622 ± 0.0003 s (where the error is given by the width 
(cr) of a Gaussian fit to the trough), in excellent agreement 
with the X-ray period of 5.7620692 ± 0.0000002 s (|Rea et al.1 
I2OO9I ). The corresponding light curves are shown as dashed 
lines in Figures [2^, ^jp and [2j; and are in good agreement 
with the solid lines obtained using the first data reduction 
technique, lending additional confidence to our reduction 
and analysis techniques. 

To test the robustness of the optical period, we 
constructed periodograms of each night separately (Fig- 
ures [3^ and[3]3); the deepest trough on each night is indeed 
coincident with the X-ray period. We also constructed 1000 
randomised light curves from the original 16 988-point light 
curve by randomly re-ordering the time-series: none of the 
resulting 1000 PDM periodograms showed a lower trough 
than the one at 5.7622 s. To check for instrumental artefacts, 
we calculated the periodograms of the sky, the comparison 
star, object position and object FWHM - none showed any 
evidence for a periodicity at 5.7622 s. As a final test, we also 
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Figure 2. Pulse profiles of SGR 0501+4516. (a) j'-band pulse 
profile on 2008 December 31. (b) i'-band pulse profile on 2009 
January 1. (c) i'-band pulse profile of both nights combined. The 
solid lines and dashed lines show the pulse profiles obtained using 
techniques (i) and (ii), respectively (see Sections l2.ll and l2.2l l. (d) 
0.5 — lOkeV pulse profile on 2008 December 24 (solid line) and 
2009 January 7 (dashed line), (e) 8-12 keV pulse profile on 2008 
September 30 (solid line) and 2008 August 29 (dashed line), (f) K- 
band pulse profile on 2009 January 28. With the exception of the 
/^-band light curve, which has had a phase ofi'set of 0.1 added (see 
Section l3.4l for details), all of the above profiles have been folded 
on the same X-ray ephemeris and hence their phases relative to 
each other are correct. For clarity, two cycles of each pulse profile 
are shown, normalised by dividing by the mean number of counts. 



calculated a sine-f itting based (Lomb-Sc argle) periodogram 
of the light curve (jPress fc Rvbickilll989l l - as expected the 
strongest peak occurs at precisely half the period determined 
from the PDM periodogram. No other significant signals at 
shorter or longer periods are present in the periodograms, 
and there is no evidence for optical bursts in the light curve 
at levels over 4(t above the mean. 



3.3 X-ray light curves 

The evidence outlined above is conclusive - we have de- 
tected optical pulsations from SGR 0501+4516 on the same 
period as the X-ray pulsations. To compare the optical and 
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Figure 3. Phase Dispersion Minimization (PDM) periodograms 
of the i'-band hght curve of SGR 0501+4516. Prom top to bottom, 
the PDM periodogram on 2008 December 31, 2009 January 1 and 
both nights combined. The extended tickmarks at a frequency of 
0.1735 Hz show the period of the X-ray pulsation. 



ULTRACAM, Smft-BAT and XMM-Newton, the accuracy 
of the absolute time-stamping of Gemini-NIRI data frames 
is not known, and hence the absolute phase of the 7f-band 
light curve is uncertain. We therefore added a phase offset 
of 0.1 (i.e. ~0.6s) to the light curve in order to align the 
two minima around phases 0.35 and 0.85 with the corre- 
sponding minima in the optical light curves. The result is 
shown in Figure [Sf. Assuming this phase offset is justified, 
it can be seen that the infrared light curve shows pulsations 
with a similar morphology to the optical light curves, al- 
beit with_a_2owerpT^^ fraction: The rms pulsed fractions 
fsee lDhillon et al.1 (|2005h for the definition used) of the light 
curves presented in Figure[2]are 52±7%, 30±4%, 26±5% and 
20±5% in the optical, 0.5-10 keV X-rays, 8-12 keV X-rays 
and infrared, respectively. 

We believe that our detection of infrared pulsations is 
robust - a periodogram analysis similar to the one we per- 
formed with the optical data (Section 13. 2|) shows a peak at 
exactly half the X-ray rotation period with a strength that 
lies in the top 0.7% of all the peaks in the periodogram. 
Although this is not significant enough to claim an inde- 
pendent detection of the period in the infrared (there are 
still many stronger peaks in the periodogram), it does lend 
support to the reality of the observed infrared pulsations. 



X-ray pulse profiles, we show 0.5 — 10 keV data obtained with 
the Swift X-ray telescope on 2008 December 24 and 2009 
January 7 in Figure Oi and 8 — 12keV data obtained with 
the XMM-Newton X-ray telescope oi l 2008 August 29 and 
2008 September 30 in Figure [2}; (see iRea et all ( |2009l ') and 
references therein). Like the optical light curves, the X-ray 
data have been referred to the Sol ar System barycen tre and 
folded on the X-ray ephemeris of iRea et al.1 (|2009l ). which 
means that the relative optical/X-ray phasing is correct. The 
0.5 — lOkeV X-ray pulse profile has a broad, single-humped 
morphology centred around phase 0.1, in stark contrast to 
the double-peaked optical profile. The broad optical peak 
occurs at approximately the same phase as the 0.5 — lOkeV 
X-ray peak, with the secondary optical maximum more or 
less coincident with the X-ray minimum. The harder X-rays 
between 8—12 keV, on other hand, look remarkably similar 
to the optical light curves. 



3.4 Infrared light curve 

As part of a study of the long-term infrared light curve 
(Levan et al., in preparation), we also obtained some high- 
speed T^-band photometry of SGR 0501+4516. These data 
were obtained with the Near Infrared Imager (NIRI) on the 
8.1m Gemini- North telescope in Hawaii on 2009 January 
28, when the SGR had a magnitu de of K = 19-7 ± 0.1 (c.f. 
K = 19.2 + 0.2 on 2008 August 25: lRea et a"Lll2008l '). The in- 
dividual frames were reduced in the standard fashion using 
the NIRI pipeline, utilizing a sky frame to remove quadrant 
offsets, and subtracting a column- stacked median in order 
to remove low-level striping in the data. Object counts were 
then extracted in the manner outlined in Section[2]2] The re- 
sulting light curve, consisting of 1168 data points each of 1 s 
exposure time and 6 s dead time, was folded into 10 phase 
bins using the same ep hemeris that was used to fold the 
optical and X-ray data (|Rea et al.1 [20091 ). However, unlike 



4 DISCUSSION AND CONCLUSIONS 

Our detection of optical pulsations on the X-ray spin 
period has confirmed that th e proposed optical an d in- 
frared counterpart r eporte d bv lFatkhuUin et al.l (|2008l ) and 
iTanvir fc VarricattI (|2008l ) is indeed SGR 0501+4516. This 
makes it the first SGR with an unambiguously-detectecQ 
counterpart at these wavelengths, providing us with a unique 
opportunity to study the emission mechanism. 

The optical properties of SGR0501+4516 are remark- 
ably similar to those exhibited by th e two AXP s with de- 
tected optical pulsat ions, 4U 0142+61 (|Kern fc Ma rtin 200i, 



Dhillon et al] |2005| ) and IE 1048.1-5937 (fPhillon et al 



2009 



All three objects exhibit optical light curves with 
broadly similar pulsed fractions and pulse profiles that are 
in phase with the X-rays. This supports the widely-accepted 
link between the AXPs and SGRs, and indicates that they 
share the same optical emission mechanism. The pulse 
profiles of SGR 0501+4516 are also remarkably similar to 
Geminga, a rotation-powered, middle-aged, radio-quiet pul- 
sar with a spin period of 0.237 s and a magnetic field strength 
approximately two orders o f magnitude lower than tha t pro- 
posed for SGR0501+4516 (|Bignami fc CaravedlT99^ '). Like 
SGR 0501+4516, Geminga exhibits double-peaked optical 
and hard X-ray pulse p rofiles and a single- peaked modu- 
lation in the soft X-rays jShearer et al.|[r998l V 

It is difficult to reconcile the results presented in Sec- 
tion [2] with the simplest form of the fall-back disc model of 
AXPs and SGRs, as any optical emissio n is predicted to be 
due to reprocessed X-rays from the disc (jPerna fc Hernguisl] 



2 SGR 1806-2 llKosugi et al.l |2005| . Ilsrael et al.1 l2005f) and 
SGR 1900+14 jTesta et al.l2008^ have faint, unconfirmed infrared 
counterparts in crowded fields, based on the observation of long- 
term variability. 
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I2000I ) and hence should have the same or lower pulsed frac- 
tion and a similar morphology. It is possible to appeal to 
beaming effects or special geometries, so that the X-ray pul- 
sations we see are different to those seen by the disc, but 
such models seem unlikely now that optical pulsations with 
similar properties have been discovered in one SGR and two 
AXPs. Further modelling would be required, however, to 
confirm/refute whether or not our observations are consis- 
tent with more elaborate forms of the fall-back disc model, 
such as the disc-star dyna mo model jCheng fc RudermanI 
ll99ll : lErtan fc Chenglliooi ). 

Pulsed optical emission in SGR 0501-f 4516 can be ex- 
plained within the context of the magnetar model by 
non-therma l emis sion from the inner magnetosphere (see 
IZane et alj (|201ll ) and references therein). In this model, 
the softest X-rays in the 0.5-12 keV range are thermal pho- 
tons from the neutron star surface, and the hardest X-rays 
in this range are produced by resonant cyclotron scatter- 
ing of the surface photons by mildly relativistic particles 
moving along the closed field lines in the inner magneto- 
sphere. The optical and infrared emission is then produced 
slightly further out in the magnetosphere by curvature ra- 
diation emitted by relativistic electrons moving along the 
closed field lines. The shared magnetospheric origin of the 
optical, infrared and harder X-rays suggest that they might 
also share similar pulse profiles, as opposed to the softer 
X-rays from the neutron star surface. This is exactly as ob- 
served in Fig ure [2] and provides strong motivation to extend 
the model of IZane et al.l (|201lf ) to predict pulse profiles for 
detailed comparison with these data. 
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